Non-local Andreev transport through an interacting quantum dot 
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We investigate sub-gap transport through a single-level quantum dot tunnel coupled to one su- 
perconducting and two normal-conducting leads. Despite the tendency of a large charging energy to 
suppress the equilibrium proximity effect, a finite Andreev current through the dot can be achieved 
in non-equilibrium situations. We propose two schemes to identify non-local Andreev transport. In 
one of them, the presence of strong Coulomb interaction leads to negative values of the non-local 
conductance as a clear signal of non-local Andreev transport. 

PACS numbers: 74.45.+c,73.23.Hk,73.63.Kv,73.21.La 



I. INTRODUCTION 

The entanglement of electrons in Cooper pairs of a su- 
perconductor can generate non-local transport effects. A 
prominent example is Crossed Andreev Reflection (CAR) 
at the contact between a superconductor and two nor- 
mal leads: there the two electrons with opposite spins 
and symmetric energies with respect to the Fermi level 
of the superconductor, that are transferred through the 
normal-superconductor interface via Andreev reflection^ 
originate from or end up in different normal leads. This 
is a non-local transport phenomenon which has been ex- 
tensively studied both theoretically^^^ 7 -^^^ and 
experimentallyii2i±£iiii±£ii& 

The main problem to identify the non-locality of CAR 
in a transport measurement is to separate it from other 
transport channels. In this paper, we study non-local An- 
dreev transport through a single-level quantum dot. The 
quantum-dot level energy can be tuned by a gate voltage, 
which opens the possibility to control the Andreev trans- 
port channels. At first glance, strong Coulomb interac- 
tion in the quantum dot seems to be counterproductive: 
the formation of a finite equilibrium superconducting pair 
amplitude is suppressed since a large charging energy pre- 
vents the equilibrium state to be a coherent superposi- 
tion of dot states with particle numbers differing by two, 
and Cooper pairs can be transferred through the dot by 
higher-order tunneling processes (cotunneling) only. On 
the other hand, a finite non- equilibrium pair amplitude 
in the dot can be achieved with a bias voltage . 17 ! 18 A 
large charging energy provides even the key ingredient 
for identifying non-local Andreev transport in one of the 
schemes we propose. 

Andreev reflection through quantum dots, a prob- 
lem which combines Coulomb interaction, superconduct- 
ing correlations and non-equilibrium, has been exten- 
sively studied theoretically i 17 d 8 i 19 i 20 i 21 i 22 i 23 i 24 Here, we 
apply the diagrammatic real-time transport theory of 
Ref. [l7d8. The relevance from the experimental point 
of view is proven by the recent success in coupling 
superconductors to quantum dots in either a carbon 
nanotub o 25 i 26 ' 27 i 28 or a semiconductor nanowire. 29,30,31 




FIG. 1: (color online) Setup: a quantum dot is tunnel coupled 
to one superconducting and two normal-conducting leads. 
The latter may be ferromagnetic with magnetization direc- 
tions riiL and riiR = ±riiL- 

In particular, we propose to investigate non-local effects 
in Andreev transport through a single-level quantum dot 
with one superconducting and two normal-conducting 
leads, which may be (i) ferromagnetic with collinear mag- 
netization or (ii) non-magnetic (see Fig. [1]). Thereby we 
identify non-locality either (i) by the dependence of the 
Andreev current on the relative orientation of the two 
ferromagnets that are biased with the same voltage, or 
(ii) by the response of the current in one normal lead to 
the applied voltage in the other one. 



II. MODEL 

The Hamiltonian of the system is H = i/ c jot+X) t; (^'? + 
#tmm,r;)- The label r/ = L,R, S corresponds to the left, 
right and the superconducting lead, respectively. The 
dot is described by the Anderson-impurity model: 

Hdot = ^ edl d a + Un^n l , (1) 

where d^ is the annihilation (creation) operator for a 
dot electron with spin a, n a = d\d a is the correspond- 
ing number operator, e the energy of the spin-degenerate 
single-particle level, and U the on-site Coulomb repul- 
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sion. The Hamiltonian of the lead 77 reads 

ka k,k' 

where the single-particle energies e^ka are spin dependent 
in the case of ferromagnetic leads (with the quantization 
axis for a being along the magnetization direction of the 
left lead), and the BCS pairing- interaction strength g v is 
non- vanishing only for 77 = S. The lead-electron opera- 
tors are c v ka and c^ fccr . We treat the superconductor on a 
mean-field level, which introduces the notion of fermionic 
quasiparticles and a Cooper-pair condensate. The cou- 
pling between the dot and the leads is taken into account 
by the tunneling Hamiltonians 

H t „ = V V Y J (clkJ* + H.c.) , (3) 

ha 

where for simplicity the tunnel matrix elements are 
considered to be spin and wavevector independent. 

The spin-resolved tunnel-coupling strengths are de- 
fined as T ncr = 2tt\V 71 \ 2 Y.k <K W _ e ka) = 2ir\V v \ 2 p va , with 
Prfa being the density of states of the spin species a in 
lead 77, which we assume to be constant. We also define 
the mean level-broadening Y ri = \ r,,^. The spin po- 
larization of lead 77 is defined as p v = {p v + — p v - )/{p-q+ + 
Pr)-), with +(— ) denoting the majority(minority) spins. 

III. METHOD 

We integrate out the leads' degrees of freedom to ob- 
tain an effective description of the quantum dot, whose 
Hilbert space is spanned by the four basis states |x) G 
{|0>, I T), I I), |D> = 44|0)}, with energies E , P T = E u 
and -Ed, corresponding to an empty, singly and dou- 
bly occupied dot. It is useful to define the detuning 
8 = Ed — E = 2e + U. The dot is described by its 
reduced density matrix, p re di whose matrix elements are 
P*^ 1 = (xi|pred|x2)- The superconducting proximity ef- 
fect induces a finite pair amplitude on the dot, expressed 
by the off-diagonal matrix element P® = (P^)* • 

In the stationary limit, the elements of the reduced 
density matrix obey the generalized master equation 

i{E Xl -E X2 )P£ = J2 W £x> P £ W 

with generalized rates W XlXl , that can be com- 
puted by means of a real-time diagrammatic 
techniquei 17 i 18 ' 32 i 33 ' 34 i 35 The stationary current out 
of lead 77 can be expressed as 

J v = -t T W XX ^P X K (5) 
" % A^i xx 2 x 2 ' w 

XX\ X2 

where W X Y' = T sW xx ? sv , and W xx } sr> is the sum of 
xx 2 L - IS XX2 XX2 

all generalized rates that describe transitions in which in 
total s electrons are removed from lead 77. 
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FIG. 2: Andreev bound-state energies as a function of the 
level position e, with Ts/U = 0.2 and Tn = 0. 

The diagrammatic rules to compute the diagrams con- 
tributing to the rates in the presence of su per conduc ting 
and ferromagnetic leads are given in Refs. Il8l and l34ll 3 5l 
respectively. In the following, we assume small and equal 
tunnel-coupling strengths IY = Tr = Tn < T to the 
left and right lead, which we keep up to first order in 
the calculation of the current. The chemical potential 
of the superconductor is chosen as the reference for the 
energies, i.e. ps — 0. To study subgap transport, we 
consider the limit of a large superconducting order pa- 
rameter A — > co, which we choose to be real. In this 
case, all orders in the tunnel-coupling strength with the 
superconductor Tg can be resummed exactly 18 i 36 In the 
absence of a superconducting lead, the excitation ener- 
gies of the dot, e and e + U, are split by the charging 
energy U. Due to the tunnel coupling to a superconduc- 
tor the particle and hole sector of the Hilbert space are 
mixed, leading to four Andreev bound-state energies, de- 
fined as poles of the retarded Green's function of the dot 
for T N = 0: 

^=4 + V( £ + ?) 2 + T- < 6 > 

with 7', 7 = ±1 (see Fig. [5]). In the following, we also 

make use of the definition £a = \J { € + -j) 2 + The 
Andreev transport is largest for small detuning \8\. We 
consider the regime \8\ < \/U 2 — FJ, for which the in- 
equality Pa,+,+ > Pa,+,- > > Pa,-,+ > Ea,-,- 
holds. 

The Andreev channel supporting transport from a nor- 
mal lead to a superconductor through a quantum dot can 
be switched between different states by an applied bias 
voltage /ZN-— In the following we consider p^ > (the 
case /in < is obtained from the symmetry transforma- 
tion /iN —> ~Pn, 8 —> —8, and Jn — * — Jn)- The different 
states are characterized by how they influence and probe 
the state of the quantum dot. 

For small bias voltage, Ea,+,- > Mn, all rates (to first 
order in Tn) involving an electron transfer from or to 
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the normal lead are not affected by the superconduc- 
tor. The current can be written in the simple form 
Jn = -tTn (1 + Q/e), independent of the pair ampli- 
tude P D , where Q/(-e) = <£ ff 4d a ) = P T + P l + 2P D 
(with — 2e < Q < 0) is the average quantum-dot charge. 
The latter can be affected by the proximity effect in 
the quantum dot. In the stationary limit, the dot is 
singly occupied, Q = — e, and the Andreev channel is 
Coulomb blocked, Jn = 0. A current can only flow 
for Q ^ — e, which, e.g., could be achieved by attach- 
ing a voltage-biased third lead. For large bias voltage, 
/iN > i?A,+,+, the Andreev channel is also independent 
of the dot pair amplitude, with Jn = Tn (1 + Q/2e). 
In both cases, the sub-gap transport involves Andreev 
processes at the interface between quantum dot and the 
superconducting lead only. As a consequence, the current 
is, in both cases, insensitive to the sign of the detuning 5. 
This is different for the regime of intermediate voltages, 
£a,+,+ > > £^A,+.-- In this case the rates involving 
an electron transfer to or from the normal lead do depend 
on S and Ts, the rates W^f and describing proxi- 
mization of the quantum dot are non-vanishing, and the 
current in the normal lead also depends on the pair am- 
plitude Pq\ An interesting feature of this regime is that 
a positive detuning drives the dot to an average occupa- 
tion of less than one electron, thus overcompensating the 
effect of the finite bias voltage^ 




FIG. 3: TMR as a function of the polarization p, for different 
values of the detuning 5 — 2e + U. Inset: Current Js in 
the superconducting lead as a function of e, for polarization 
p — 0.5 and parallel and antiparallel alignment of the leads' 
magnetizations. For both graphs, we chose /il = /Ur = Hn = 
U, /i S = 0, Ts = 0.2U, r L = r R = 0.00H7, and k B T = 0.01(7. 



of p the TMR exhibits a quadratic dependence on p, 

1 - n 7=± (E y =± 77n(£a, y , 7 ))] r 2 P 2 



TMR: 



4f 2 



E 7 ,y=± 77'(e + f + e A )f N (E A ^ <1 ) 



IV. CROSSED ANDREEV REFLECTION 



which simplifies in the regime of large bias voltage to 
TMR « p 2 for \S\ < Tg and TMR w p 2 T 2 /8 2 for \S\ > 
Ts, and in the regime of intermediate (positive) bias volt- 
age to TMR « 4p 2 /3 for |5| < T s , TMR « p 2 T 2 /S 2 for 
-S > T s , and TMR « 2p 2 for 5 > T s . 



We first consider the case of ferromagnetic leads 
with equal polarization strengths \pi,\ — \pn\ = p, 
kept at the same chemical potential /iN and char- 
acterized by the same Fermi distribution /n(c) 
1 + exp (e — /xn)Ab T]~ ■ Crossed Andreev transport is 
identified by its dependence on the relative orientation 
of the ferromagnets, quantified by the Tunneling Mag- 
neto Resistance (TMR) TMR = (J^ p - Jg)/J£ p , where 
jp~ r > = (2er s /7i)ImP^ is the current in the supercon- 
ductor for parallel (antiparallel) alignment of the mag- 
netizations. The pair amplitude P® — (e?x^t) ^ s ca l cu ~ 
lated in the stationary limit from the master equation 

Po = E x ^0x Xp x / [ w oa° - *( 2e + u )] ■ In the inset of 
Fig. [3]we show the current in the superconducting lead for 

both the parallel and antiparallel alignment of the mag- 
netizations as a function of the gate voltage for /iN = U 
(large-bias regime). The current shows a peak around 
zero detuning, e = —U/2, with a width given by r$. 
The sub-gap current for the parallel alignment is clearly 
suppressed, indicating the presence of CAR. The TMR is 
plotted in Fig. [3] as a function of the polarization strength 
p for different values of the detuning 5. For small values 



V. NEGATIVE NON-LOCAL CONDUCTANCE 

We now look for non-local effects in Andreev transport 
for non-magnetic leads, pl = PR = 0, by studying the 
current in the right lead, Jr , as a function of the voltage 
applied to the left lead /xl- For a three-terminal device, 
we define a non-local conductance, G nl = Jr/(^s — Ml), 
as the current response in the right lead to a voltage 
bias between superconductor and left lead. In particu- 
lar, we will consider the non-local differential conduc- 
tance GJjjff = —dJji/dfij^. Direct transport between 
the two normal leads contributes with a positive sign to 
G^] ff . Non-local transport channels such as CAR may 
contribute with a negative sign. For a large class of such 
three-terminal devices, however, it has been shown^i that 
the sum of all contributions to G^\ s remains positive. In 
contrast, we find for our system regimes with negative 
values of GJjjg and, even more striking, negative values 
ofG nl . 

In Fig. 3^,, b, and c, we show the current in the right 
lead as a function of the electrochemical potential in the 
left lead for the three cases of a small, intermediate and 
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FIG. 4: (color online) Current Jr in the right lead as a func- 
tion of the chemical potential [il of the left lead for different 
values of the detuning 8 and T s = 0.2C7, IY = Tr = O.OOlt/, 
Pl = PR = 0, ksT = O.Olf/, and (a) fin = (small bias), (b) 
fi>n = 0.5(7 (intermediate bias), and (c) fj, R = U (large bias). 



Figs. |4^ and c, the plateau height of the current, given by 
J R = 2er N (1 + Q/e) and Jr = 2er N (1 + Q/2e), respec- 
tively, directly reflects the average quantum-dot charge, 
which is strongly influenced by the proximity effect i 18 ! 38 
For (5 = 0, the proximity effect is maximal, with Q = —e 
for all values of the bias voltages, which leads to Jr — 
and Jr = eF^ in Fig. QJi and c, respectively. With in- 
creasing \S\, the proximity effect decreases, and the cur- 
rent approaches the value expected in the absence of the 
superconducting lead. Third, the striking feature indi- 
cating non-local Andreev transport is the non-monotonic 
dependence of Jr on /xl, i.e., the appearance of a nega- 
tive non-local differential conductance, G^\ s < 0. Even 
more remarkable is for //r = the negative non-local 
conductance G nl < that occurs for positive/negative 
detuning S at positive/negative in the intermediate- 
bias regime. To understand this behavior, we realize 
that in this regime there are combined Andreev processes 
that involve both the interfaces from the quantum dot 
to the superconductor and the left lead, while there are 
no Andreev processes involving electron transfer between 
quantum dot and right lead. The intermediate-voltage 
Andreev transport between left lead and superconductor 
yields an average dot charge that is determined by the 
dot level's position relative to the chemical potential of 
the superconductor rather than that of the normal lead: 
for positive (negative) detuning the probability of double 
occupation decreases (increases), and the average occu- 
pation of the dot is smaller (larger) than one. This devi- 
ation from single occupancy is probed by the right lead. 
Changing the sign of detuning leads to a sign change in 
the current measured in the right lead. 

We remark that the negative non-local conductance is 
not due to CAR. In fact, making the normal-conducting 
leads ferromagnetic suppresses the negative non-local 
conductance, independent of whether the ferromagnets 
are aligned parallel or antiparallel, since a finite spin ac- 
cumulation on the quantum dot reduces the dot pair am- 
plitude. For CAR an enhanced (reduced) effect would be 
expected for the antiparallel (parallel) alignment. The ef- 
fect that we predict is rather a consequence of combined 
Andreev processes between left lead and superconductor. 
The negative non-local conductance can only be probed 
because of a large charging energy that prohibits direct 
transport between the normal-conducting leads. 



VI. CONCLUSIONS 



large voltage /xr applied to the current probe, respec- 
tively. Several features of these current-voltage charac- 
teristics indicate the strong coupling of the quantum dot 
to a superconducting lead. First, there are four steps (in- 
stead of two) associated with the four Andreev bound- 
state energies. Second, the height of many of the plateaus 
is sensitive to the detuning 6. In the cases displayed in 



We investigated non-local Andreev transport through 
an interacting quantum dot in a three-terminal setup 
with one superconducting and two normal-conducting 
leads. We considered two different biasing schemes. In 
the first one, the normal-conducting leads are ferromag- 
netic with collinear magnetizations and they are kept at 
the same chemical potential. The key results for this case 
is that CAR occurs due to the non-equilibrium proxim- 
ity effect in the dot and it is characterized by a finite 
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TMR. In the second scheme, the two normal-conducting 
leads are non magnetic, and the response of the current 
in one lead to the voltage applied to the other one is 
studied. In that case, non-local Andreev transport is 
identified by negative values of the non-local differential 
conductance. Even more strikingly, we find regimes with 
negative values of the full non-local conductance. The 
virtue of employing quantum dots lies, first, in the pos- 
sibility to tunc the Andreev channels by a gate voltage 
and, second, in the presence of a large charging energy 
which generates specific transport regimes characterized 



by a negative non-local conductance. Both aspects are 
advantageous for a clear identification of non-local An- 
dreev transport. 
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